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Abstract 

Experimental measurement of the elliptic flow parameter V2 and hydrodynamic model together 
showed that thermalization in the central region at the Relativistic Heavy Ion Collider to be 
perplexingly fast. This is a mystery in itself since none of the numerical perturbative QCD models 
are able to achieve such a feat. By exploiting a theoretical oversight on collinear processes in an 
out-of-equilibrium system it is argued that, in the bottom-up thermalization scenario, equilibration 
can proceed at a higher rate than what is expected in the conventional perturbative QCD picture. 

PACS numbers: 25.75.-q, 12.38.Bx, 12.38.Mh, 24.85. +p 
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I. INTRODUCTION 

Recently measurement by STAR QQQ, PHENIX Q and PHOBOS collaborations 
of the elliptic flow parameter f 2 in the central region at the Relativistic Heavy Ion Collider 
(RHIC)at Brookhaven agreed with that generated from simulation using the hydrodynamic 
model ^7,0,0]. As is well-known that hydrodynamical model required that a system to 
be in complete equilibrium. This implies that the central region in the collisions at RHIC 
achieved equilibrium within an extremely short time of around 0.6 fm/c. In case one has 
any doubts about this conclusion, the authors investigated an alternate scenario where the 
collision system was thermalized only in the transverse directions We flnd that it is 
impossible to achieve the same V2 results while keeping the same flnal transverse momentum 
distribution. Even trying to maintain the latter required very drastic modification to the 
initial conditions so that these become unrealistic. This is both good and bad news from 
a theoretical point of view. It is good because a fully thermalized environment at least in 
the central region permits a lot of simplifications in the calculation. A thermalized system 
is much simpler than an out-of-equilibrium system. However this is also bad news because 
it shows that we have not fully understood what really happens in the initial stage of the 
collisions. There are two reasons for that. First it was shown in 111, 12, 13| that the special 
case of the pano„ cascade .ode. Q H Q w.h o.lv e.a.,c collision. .e,u.ed either 
an initial parton density or the elastic cross section 15 times greater than the actual value 
in order to reproduce the same elliptic fiow measured in the experiments. This is a huge 
factor, one does not expect, based on today's knowledge, that including inelastic collisions 
will produce a large enough total cross section. Second it is our opinion that none of the 
parton based numerical models for simulating heavy ion collisions, parton cascade included, 
is showing such perplexingly fast thermalization |l4 . Iisl IigI Q| . They do show signs that 
the system is approaching equilibrium at a rapid pace even in the face of the longitudinal 
expansion, but none of these models managed to thermalize within as short an interval 
as 0.6 fm/c. We do not consider the results from other non-parton based models because 
they use other less well understood and less rigorous mechanisms such as strings, hadronic 
collisions etc or the combinations any of these. In any case it is hard to justify the use of any 
hadronic mechanisms in view of the energy density that can be reached at RHIC energies. 
There is one other mechanism that can be at work which is that as the system cools, the 
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average energy of the system must come down. This will lead to an progressive increase in 
the strength of the interactions This mechanism can certainly help with equilibration. 
However this mechanism does not come into play until at a much later stage in the collisions, 
well beyond the 0.6 fm/c limit imposed by the V2 result. 

To solve this very rapid equilibration problem, the first question that one has to decide 
is what mechanism leads to this fantastic result. Is it perturbative or nonperturbative in 
nature? Parton cascade model is supposed to have incorporated all known and relevant 
perturbative QCD effects One can easily conclude that the as-yet- unknown mechanism 
must be nonperturbative. Nevertheless another picture of the nuclei exists that involved 
viewing the initial nuclei as a frozen crystalline glass of QCD color Initial conditions 
produced from this model have been used to investigate and tried to answer the perturbative 
or nonperturbative equilibration question This line of investigation resulted eventually 



in ref. 



2l| , this paper used a result already pointed out in ^ which is that thermalization 



cannot come from leading order processes alone, higher order processes are important too. 



It was concluded in 



2l| that the equilibration could proceed perturbatively. However this 



paper did not address how rapid a perturbative QCD based equilibration could proceed. 
Even if equilibration in heavy ion collisions can be achieved via perturbative QCD alone, 
can it be done with such phenomenal speed as revealed by the combination of elliptic flow 
measurements and the hydrodynamic model? There is a clear advantage of a perturbation 
theory based rapid equilibration as opposed to a nonperturbative one. The former is well 
understood theoretically and tested in experiments while the latter is less rigorous and based 
on a much less solid foundation. Therefore before one draws any definitive conclusion based 
on the results of the parton based models and gives up on perturbative QCD, another careful 
examination is warranted. 

We will show in this paper that there is indeed a new mechanism that has been overlooked 
largely because of a piece of well-known knowledge that has been overly generalized and 
assumed to hold under all circumstances even though no careful study has been done in 
some situations. By this we means the cancellation of coUinear divergence for a theory 
of massless particles which has been shown to be true in the vacuum 23|, |2J] and in an 
equilibrium system. Indeed this is true in the vacuum as total cross-sections are free of any 
coUinear logarithms as these are cancelled between the real and virtual contributions. This 
is also true for a system of massless particles in equilibrium. Because of this, one tends 
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to assume tacitly that the same will also hold in a non-equilibrium environment. This is 
not true. In general this cancellation between real and virtual graphs ceases to occur for 
a system that is not in equilibrium. One might be alarmed by this statement since there 
should not be any remaining divergences. This is reasonable but thanks to screening in 
a medium such divergences do not exist. The cancellation of collinear logarithms for a 
system in equilibrium is thus a double safeguard against this kind of divergences. In such a 
system, removing either screening or the cancellation of collinear logarithms alone does not 
result in any divergences. In the rest of the paper, we will show in details how there are 
necessarily non-cancelled collinear logarithms in an out-of-equilibrium parton plasma and 
how this result can be exploited to bring about a thermalization via perturbative QCD that 
can proceed at a higher speed than any parton models have shown thus far. 

II. THE MOST RELEVANT COLLISION PROCESSES FOR EQUILIBRATION 

In the early stage of the system, the plasma is gluon dominated because of the small-x 
growth in the gluon distribution in the nucleons and also because of the stronger gluon- 
gluon interactions than that of gluon- (anti) quark or (anti) quark- (anti) quark. Therefore for 
our present attempt at explaining the very fast thermalization seen in experiments, one only 
has to consider a pure gluon plasma. In general there are many possible interactions even 
only among gluons. Not all of them are useful for the purpose of equilibration. We will 
group the different contributions below and discuss the importance and relevance of them 
with regard to thermalization. 

A. 2 < — > 2 scattering 

(i) Small angle collision: At the leading order, it is the simple 2-to-2 gluon-gluon scatter- 
ing. It is well known that this process is dominated by small angle scatterings because 
of the familiar Coulomb exchange divergence at small momentum transfer. Although 
the probability of this process is large, the outgoing momenta are not that different 
from the incoming ones since they are only rotated somewhat in momentum space 
from the incoming momenta. Thus this normally dominant process is in fact terribly 
inefficient at driving the system towards equilibration. 
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(ii) Large angle collision: At the same order, there is the remaining large angle scattering 
process. The outgoing momenta are necessarily very different from the incoming ones 
therefore large angle scatterings are much better at redistributing momenta which is 
the essence of what thermalization is about. The drawback of this process is that 
the probability is not nearly as large as the small angle collisions. There will be an 
occasional large angle collision amongst the much more frequent small angle collisions. 

At this level of the basic 2-to-2 scattering, the possibility for the driving force behind the 
equilibration process is restricted between a probable but inefficient small angle scatterings 
and a less probable but more efficient large angle collisions. This is not very promising. 
Indeed as have been pointed out in the second reference of the elastic scattering process 
is not the main driving force for equilibration. This is because in that numerical study, the 
entropy generation rate from 2-to-2 scatterings is not only comparable but sometimes even 
subdominant to that of the 2-to-3 gluon processes. This point of view is reinforced in a later 
paper [2^. Therefore higher order processes must be considered. The key question is which 
higher processes? 

B. 2 < — > 3 scattering 

We now attempt to identify what processes at this order are most relevant and most 
important for thermalization. 

(i) Small angle collisions with a small angle gluon radiation/absorption: At order there 
is the possibility of a gluon emission on top of the main scattering. Since the emitted 
gluon can be of any magnitude limited only by energy-momentum conservation, this is 
a much more flexible process for achieving equilibration. However flexibility is not the 
only factor that one has to consider. In this case, all the momenta are aligned around 
the directions of the two incoming momenta. So there is still no significant momentum 
transfer here. Besides there is the extra power of as that reduces the effectiveness of 
this contribution. The same applies to the small angle absorption. 

(ii) Small angle collisions with a large angle gluon radiation/absorption: This type of 
contributions come in two different forms. One can be reclassified as large angle 
collisions with a small angle emission or absorption to be discussed below. Those 
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contributions that cannot be reclassified are as before down by a power of as although 
in view of the large angle gluon emission and the small angle scattering, they are 
potentially useful for achieving equilibration. 

(iii) Large angle collisions with a small angle gluon radiation/absorption: Again large angle 
collisions are important for sizable momentum rearrangement as already mentioned, 
what is more important is the well-known fact that small angle gluon emission or 
absorption accompanying a hard collision can give rise to a large collinear logarithm 
in the probability provided that the collision is sufScently hard. In the medium this is 
determined by both the hardness of the collisions and the medium screening, this will 
be seen in later sections. These kind of large logarithms can compensate for the small 
power of as that comes with the gluon emission or absorption. Thus the contribution 
from this type of processes is comparable in size to that from the leading order. 

(iv) Large angle collisions with a large angle gluon radiation/absorption: Large angle pro- 
cesses sadly do not receive any collinear logarithmic compensation of the small cou- 
pling. Processes discussed in (iii) are definitely more important for this reason. 

C. Higher order processes 

Going to even higher orders, one can have various mixtures of large and small angle 
processes given that the possibilities are nearly endless as one goes to higher and higher order. 
After all we are using perturbative QCD, higher processes are in general less important. 
Prom the above discussion, we have ascertained that large angle processes are useful for 
equilibration because of their ability to rearrange the momentum distribution. However they 
are weighed down by the coupling and have no singularity to enhance the frequency of their 
occurrence. It is certain therefore that there cannot be too many large angle subprocesses 
in any given scattering and there cannot be none either. The simplest way to proceed is to 
combine both large and small angle processes. Since there must be at least one large angle 
process and the coupling associated with small angle subprocesses can only be compensated 
by a large collinear logarithm if there is a hard scattering and/or small screening, it is logical 
that for contribution at any order to have one hard collision to be accompanied by all sorts 
of collinear subprocesses. These type of contributions are now comparable in size to that of 
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the leading order 2-to-2 scattering. These are the processes that we considered to be most 
important for equihbration and we will concentrate on these contributions and consider 
them further up to order in the rest of this paper. Note that we are not relying solely on 
the momentum rearranging power of the large angle collisions. They serves mainly to direct 
and guide where the remaining gluon radiations will populate the momentum space. It is 
these radiations that actually help to bring about thermalization of the system.. In a way we 
are arguing for usirig multi-gluon processes to achieve rapid thermalization not dissimilar to 
what was done in j^l for achieving chemical equilibration. Whereas in [2^ they included 
mostly every ng ^ mg process indiscriminately, we consider only what in our opinion the 
most important processes for thermalization here. 

Notwithstanding of what was mentioned in the introduction, any alert readers will rightly 
object that one cannot use coUinear processes in this way because of the cancellation between 



the real and virtual contribution 



23[. This is true only if there is a sum over degenerate 



states j^]. The equilibration of a gluon plasma, and in fact for other many-body systems as 
well, is not at all inclusive because momentum regions are separate from one another. The 
processes cannot be any more exclusive than this but it is usual to study interactions and 
momentum states collectively, in which case one has indeed to consider this cancellation. As 
mentioned in the introduction and we will see below, this cancellation does not hold in an 
ou.of-e,uiHb..„ syste. Q. But the. is „o „eed to wo.,, about any divergence because 
where there is screening, the interaction probability is always finite. 



III. THE LEADING AND NEXT-TO-LEADING ORDER COLLISION PRO- 
CESSES 

An out of equilibrium system is necessarily time dependent, there must be a transport 
equation to study such a system. The simplest and possibly the most used is the Bolzmann 
equation 

P^^^ = ^(-'P)' (1) 
On the right hand side are the collision terms which contain all the microscopic interactions. 
Since they are the main reason that a many-body system can equilibrate, we will spend the 
rest of the article studying these. As is commonplace in perturbative QCD, the collinear 
gluons are either on-shell or are nearly on-shell so that only the transverse degrees of freedom 
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FIG. 1: The basic scattering of gluon 1 with gluon 2 resulting in gluon 3 and 4. Likewise one can 
view this equahy as gluon 3 scatters with gluon 4 giving gluon 1 and 2. 

are relevant. Therefore a physical gauge such as the lightcone gauge will be used throughout. 



A. The leading gluon-gluon collision process 

The basic gluon-gluon collision is well-known. Graphically we represent it by Fig. ^ The 
collision term is 



12 



112->34| 

X[/l/2(l+/3)(l+/4) - (l+/l)(l+/2)/3/4] . (2) 

Here fj = f{pj), [dp] represents the invariant momentum-space integration measure for an 
on-shell gluon with momentum p, 

(3) 

with Ep = \p\ for massless gluons, and -Di2,34 is the energy-momentum conserving (5-function 

^12,34 = (27r)^5 W (pi +p,-p,-p,). (4) 

The square of the matrix element is given by [2^ 

^Ml^ = ^ E l-^--34r = 72 u^aliQ^) (3 - ^ - ^ - J) (5) 

where the sum is over spins and colors, s = (^1+^2)^! t = {pi—psY, and u = {pi—piY = — s — 
t. For the scale of the running coupling constant ^^(Q^) we can take = min(|s|, \t\, \u\). 
At momentum scales where there are three active quark flavors, it can be parametrized as 

47r 

/?oln(Q7L(3)) 

where L(3) = 0.246 GeV and (3o = 9. We will assume that Q is significantly larger than the 
scale L associated with non-perturbative QCD effects. 
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B. The gluon-gluon collision with coUinear emission or absorption 



For 5-gluon processes there are two distinct contributions to the colhsion terms: 

x[/l/2(l+/3)(l+/4)(l + /5) - (l + /l)(l + /2)/3/4/5], (7) 

and 

Cl23^45(Pl) = y C?$23C?$45^123,45 l-M 123^451 

x[/l/2/3(l + /4)(l + /5) - (l + /l)(l + /2)(l + /3)/4/5]. (8) 

We have introduced a compact notation for the multiparticle phase space: 

J d^i2...m = [dpi] [dp2] ■ ■ ■ [dp„,] . (9) 

Every such factor is accompanied by a symmetry factor 1/m!. The energy-momentum 
conserving delta functions -Di2,345 and -Di23,45 are the obvious generahzations of Eq. (jH). 

Suppose the 5-gIuon processes 12 ^ 345 in Eq. (|7j) and 123 <-> 45 in Eq. (jSJ consist of a 
hard-scattering collision involving four gluons. Then two of the incoming or outgoing gluons 
can be hard but collinear. If the two gluons are both in the initial or both in the final 
state, collinear means that their invariant mass Sij = {pi +Pj)'^ is significantly less than Q^, 
where Q is the momentum scale of the hard scattering. If one is in the initial state and 
one is in the final state, collinear means that their momentum transfer \tij\ = —{pi —PjY is 
likewise significantly less than . For the process 12 <-> 345, there are nine possibilities for 
the collinear pair of gluons By exploiting the Bose symmetry, we can take the collinear 
pair to be (1,5), (2,5) or (4,5). For the process 123 45, there are again nine possibilities 
for the collinear pair. Again by exploiting the Bose symmetry, we can take the collinear pair 
to be (1,3), (2,3) or (4,3). For each of the possible collinear pairs, we wish to calculate the 
contribution to the collision terms C{p) from the collinear region of phase space. We will 
label the contribution in which gluons i and j are collinear by i || j. 

1. Both collinear gluons in the initial or final state 

We first consider the case in which the collinear gluons are either both in the initial state 
or else both in the final state and neither is the external gluon 1 (of course there is nothing 
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special about gluon 1 except its momentum has been chosen to be the argument of the 
coUision terms). We begin by considering the scattering process 12 ^ 3ij in which the two 
colhnear gluons are both in the final state with momenta Pi and pj. The matrix element 
M i2-*3ij has a pole in the invariant mass Sij. In the collinear region in which Sij <^ Q , we 
can interpret the process as proceeding through a hard-scattering 12 34* that creates a 
virtual gluon 4* with momentum p^* = Pi + pj followed by the decay of the virtual gluon 
into the two collinear gluons i and j. 

It is convenient to introduce light-cone variables defined by the direction of the 3- 
momentum p^* = n + pj in the rest frame of the plasma. The momenta of the collinear 
gluons and the virtual gluons can be written as 

Pi = (zpt,^,p±) (10) 
V zpj / 

P. - ((1-.).|,^^^,-P.) (11) 

P4. = (pt,'^,oA (12) 

\ pj y 

where p'l is the total light-cone momentum of the collinear gluons. The variable z is the 
fraction of that light-cone momentum that comes from gluon i. The light-cone energies are 
determined by the mass-shell constraints pf = p'j = 0. The invariant mass Sij of the virtual 
gluon is determined by conservation of light-cone energy: 

2 2 2 

% = ^ = -T, N ■ (1^) 

z 1 - z z{l - z) 

The phase space integral in light-cone variables is 



dp~^ d?p±_ dz d'^p± 
A7TP+ (27r)2 ~ Attz (27r) 



up u p± u,^ u, p± 

[dp] = TTZTTT^ = — ' (14) 



the latter form applies when p is the momentum of the daughter of a parent gluon from which 
a fraction z of its light-cone momentum has been inherited. We define the momentum p^ of 
an on-shell gluon that has the same light-cone momentum as the virtual gluon 4* but whose 
light-cone energy is determined by the mass shell constraint p| = 0: 

P4= {pt,0,O±) . (15) 

We proceed to make approximations to the collision integral Eq. (|7j) (with 45 replaced by 
ij) that are valid in the collinear region. In the amplitude for 12 34*, we can replace the 
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FIG. 2: An example of the 2-to-3 or 3-to-2 gluon process when gluon 4||5. 

momentum p^* of the virtual gluon by the quasi on-shell momentum p^, since the invariant 
mass Sij is small compared to the hard-scattering scale. In the amplitude for the decay 
4* —>■ ij, we keep only the most singular term which will give rise to a logarithm of Q after 
integrating over phase space. The square of the matrix element, averaged over initial spins 
and colors and summed over those of the final state, then reduces to 

E |A<12-.3./ - E |-M,,^34r (16) 

Sij 

where Pgg{z) is the unregularized g ^ gg splitting function 

where Ca = 3 is a color factor. The scale of the explicit factor of as in Eq. ()lfi|l is y/Sij- We 
can make similar approximations in the phase space integrals. We first express the phase 
space for gluons i and j in an iterated form involving an integral over the momentum of the 
virtual gluon 4* 

In the collinear region, we can approximate p^* by p^ in the measure [dp^*] and in the delta 
function -Di2^34* ■ The expression for the phase space integral Eq. (fTHj) then reduces to 

j [dpi\[dpj\D 12,^, j ~ j [dpi]Di2,u ^Q^2^^^t^ = J MP4] ^12,34 ^Y^J^ • (19) 

Combining the approximations Eqs. (fTBj) and (fT^ and using Eq. (fT^. we get 

I rf<f3.,.Dl2,3.,$]|A^12^3..f ^ I rf$34/^12,34$^|A^12^34pX j ^^^^P,,{z). 

(20) 

Our final step is to use a collinear approximation in the distribution functions fi and 
fj for the gluons. Since their momenta have components of order Q, we can neglect their 
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transverse momenta and approximate them by f^iz) and /4(1 — z), where 

f,{z) ^ f{zp,) . (21) 

Thus our final expression for the contribution to the colhsion integral Eq. ((7j) from the region 
in which gluons 4 and 5 are collinear is 

Ug J J S45 Ztt 

X /l/2(l + /3)(1 + /4(^))(1 + /4(1 - Z)) - (1 + + f2)hh{z)U{l - Z) . 

(22) 

Depicted in Fig. Elis the feynman graph for this contribution when 4||5. The contribution 
from the collinear regions 3 || 4 and 3 || 5 are given by exactly the same expression after 
appropriate relabelling of the momenta. 

The contribution to the collision integral Eq. (jSJ from the collinear region 2||3 is given by 
a similar expression. After relabelling the dummy momentum variables 45 — ^ 34, it can be 
written as 



Cll3^45(Pl) - /Mrf$34A2,34^$:i-^12^34rX | ^^f^^^.l^) 
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hf2{z)h{\ - Z)(l + /3)(1 + /4) - (1 + /l)(l + h{z)){\ + /2(1 - ^))/3/4 

(23) 



In both Eqs. ()22j) and (jzHl) . the ranges of the integrals over z and Sij have not been specified. 
The ranges in which our collinear approximation remain valid will be specified later after a 
discussion of hard-thermal loop corrections. 



2. One collinear gluon in the initial and the other in the final state 

We next consider the case in which one of the collinear gluons is in the initial state 
and one is in the final state and neither is the external gluon 1. We begin by considering 
the scattering process li 34j, where gluons i and j are collinear. The matrix element 
-^ii^34j has a pole in the variable tij = {pi — pj)"^. In the collinear region \tij\ ^ Q^, we 
can interpret the process as the splitting of gluon i into a real gluon j and a virtual gluon 
2* with momentum p2* = Pi — Pj, followed by the hard scattering 12* —>■ 34. 
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FIG. 3: An example of the 2-to-3 or 3-to-2 gluon process when gluon 2||5. 

We introduce light-cone variables defined by the direction of pi in the rest frame of the 
plasma. The momenta of the coUinear gluons and the virtual gluon 2* can be written as 

P^ = (24) 
f {l-x)pt xpI \ 

P2* = (pt,'-'^,P±) (26) 

where is the difference between the light-cone momenta of the coUinear gluons. The 
variable x is the fraction of the light-cone momentum of gluon i that is carried by the virtual 
gluon 2*. The light-cone energies of i and j are determined by the mass shell constraints 
Pi = p] = 0- The invariant mass of the virtual gluon is determined by the conservation of 
light-cone energy 



(27) 

We define the momentum p2 of an on-shell gluon that has the same light-cone momentum 
as the virtual gluon 2* and the same transverse momentum as gluon i: 

P2 = (p^,0,0x). (28) 

We proceed to make approximations to the collision integral Eq. ((7j) (with the pair (2,5) 
replaced by (ij)) that are valid in the coUinear region Fig. El The square of the matrix 
element can be approximated as follows 



E - ^^^^^ E l-^--4r (29) 

where the scale of the explicit factor of is ^f—Uj- The phase space integral can be 
approximated by 

j m[dp,]D,,,,, ^ j [dP2]Dr.,u J^2^f^^^ = I [dp2]D,2,,,^^^tpl . (30) 
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The distribution functions and fj can be approximated by /2(^) and /2(^^), respectively. 
Combining all these approximations, the final expression for the contribution to the collision 
integral Eq. (0) from the region 2||5 is 

X /l/2(^)(l + /3)(1 + /4)(1 + /2(^)) - (1 + /l)(l + /2(^))/3/4/2(i^)" • 

(31) 

The identical expression is obtained for 2||3 or 2||4. 

We can derive a similar expression for the contribution to the collision integral Eq. (jH} 
from regions in which one of the coUinear gluon is in the initial state, the other is in the 
final state, and neither one is pi. The contribution from the region 3||5 can be written as 
after relabelling of momenta 

C!!LM - /Wrf$34A2,34^El-^l^-34pX / dx^^^P,,{x) 

/l/2(l + h)f^{^){l + /4(^)) - (1 + /l)(l + /2)/3(l + /4(^))/4(^) 



X 



(32) 

An identical expression is obtained, after relabelling, for the regions 2||4, 2||5, and 3||4. In 
Eqs. ()31|) and (jH^ . the ranges of x and tjj in which our collinear approximations remain 
valid will be specified later after a discussion of hard-thermal-loop corrections. 



C. Hard-thermal loop corrections 

One important class of higher order corrections to the 4-gluon collision integral Ci2-+34(pi) 
is hard-thermal-loop (HTL) corrections to the propagators of the four hard scattered gluons 

nn 

|27L I28| . These corrections arise from hard gluons in the medium that interact with one of 
the hard-scattered gluons but are then scattered back into their original momentum state. 
These forward scattering contributions must be added coherently over the distribution of 
hard gluons in the medium. The effect of these HTL corrections on transverse gluons with 
momenta near the light cone is very simple. They become quasi-particles with a mass mt 
given by 

dppf{p) (33) 



24a, (/i^ 



TT 
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where f{p) is the distribution of hard gluons in the medium. An appropriate choice for the 
scale fi of as is the momentum scale that dominates the integral. As an estimate for the 
scale, we will use 

^2 _ irjpp''f(p) _ (34) 

In the case of the equilibrium Bose-Einstein distribution, this description gives /i = 
We will usually refer to the transverse gluon quasi-particle simply as gluons and to the mass 
given by Eq. plj) as the transverse mass. The transverse mass is related to the Debye 
screening mass mn simply by rrit = mo/V^ (sQ^. This relation arises from the structure of 
HTL corrections to the gluon propagator which are strongly constrained by gauge invariance. 

The HTL corrections will have particularly important effects in the regions of phase space 
where gluons become collinear. In particular, they provide infrared cutoff on the integrals 
over the momentum fraction and invariant masses that appear in the collinear contributions 
to the collision integrals. We will use the HTL corrections to estimate the boundaries of the 
integration regions where our collinear approximations remain valid. 

We first consider the integral over z and S45 in the 4 || 5 contribution to the collision 
integral 6*12^345 which is given in Eq. ()22|). One effect of the HTL corrections is to replace 
the propagator 1/s of the virtual gluon by — m^). Another effect is to change the 
mass-shell constraints on the collinear gluons to pi = pi = mf. This changes the expression 
Eq. (jl3p for the invariant mass of the virtual gluon to 

For fixed z, the lower limit on S45 occurs at p± = 0: 



(36) 



111-111 / 1 \ - 

z{l - z) 

The upper limit on s is set by the scale of the hard scattering. We choose to write it as 

Smax = (?Q^ + (37) 

where c is a constant near 1. The change from varying c by a factor of 2 should be included 
in the theoretical error. We will take the virtuality S45 — of the virtual gluon to also 
provide the scale of the factor as in the collinear correction factors. The integral over S45 
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can then be evaluated analytically 



ds a.js - mX) _ 2 ^ / 1^ {fQ^jl}] ^ ^^^^ 



The condition Smax > Smin constrains the integration region for z to the range z^ax > z > 



z^in, where 



1 (l + v/l-4m2/(c2g2 + ^2)j (39) 
I (l - - 4m2/(c2g2 + ^2) j _ (40) 



The condition ^^ax > ^min also imposes a lower bound on Q 

12 ^ 2 _2 



2. 



> . (41) 

We next consider the integral over x and t25 in the 2 || 5 contribution to the collision 
integral Ci2->345, which is given in Eq. ()3ip. The HTL corrections replace the propagator 
1/^25 of the the virtual gluon by 1/(^25 "~ ^t)- They also change the mass shell constraint to 
P2 — Pt — The resulting expression Eq. (f?fj) for the invariant mass of the virtual gluon 
is replaced by 

= . (42) 

1 — X 

For fixed x, the lower limit on —^25 occurs at p± = 0: 

2 

^min Z ) (43) 

I — X 

the upper limit on —^45 is set by the scale Q of the hard scattering. We choose to write it as 

- tmax = c'g' - ml (44) 

where c is again a constant near 1. The change from varying c by a factor of 2 should be 
included in the theoretical error. We will take the virtuality —^25 + mf of the virtual gluon 
to also be the scale of as in the coUinear correction factor. The integral over t25 can then 
be evaluated analytically 

d{-t) as{-t + m?) _ 2 J / (^'^'/^'. , _ (45) 



-t + mt 271 Po I ln((2-x)m2/(l-x)L 
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The condition — tmax > — ^min constrains the integration range for x to the region x < Xmax, 
where 

2 

a:max = 1 - ^ • (46) 

For the lower hmit, we will set 

The condition Xmax > imposes a lower bound on Q'^ 

1 



> . (48) 



Note that in the limit where Q'^/m'^ ^ 1, there is little distinction between the different 

limits since 2^max — ■^max and Z^am — 3^min- 



D. The collinear virtual contributions at the next-to-leading order 

We have not considered virtual corrections to the basic hard scattering discussed in Sec. 
mil These are of the same order in as as the 5-gluon processes. In the collinear limit, 
these corrections can be written in a compact way. To see this, recall that in the vacuum 
theory if one starts out with a virtual gluon, it has a chance of splitting into two gluons with 
momentum fractions z and 1 — z. The associated variation in the probability density is [s^ 

V,, + dPgg = 5(1 -z) + ^P;,iz)ds (49) 

while the gluon reduces the square of its invariant mass from s to a value between s and 
s — ds for a timelike branching. For spacelike branching, one substitutes dt for ds and t 
is changed to a value between t and t + dt. P^g{z) is the usual kernel of the Dokshitzer- 
Gribov-Lipatov-Altarelli-Parisi (DGLAP) equation 3Cl 

function part. It is the 'plus' distribution regularized gluon splitting function 

momentum is always conserved, the total momentum must always be equal to that of the 

original virtual gluon independent of whether splitting occurred, then P^g{z) must satisfy 

the identity (in the absence of q and q) 

-1 

dzzP^g{z) = (50) 



or if we write out Pgg{z) separately the parts with and without the 5{1 — z) function 

Pa,{z) = PMW + Ct=o 5(1 - z) , (51) 
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32] and it carries also a 5(1 — z) 



3G 



Since 



where Ct=o is the constant in front of the ^-function part of P^g{z), momentum conservation 
gives 

^ dz\^zPgg{z)\+ + Ct=o 6(1 - ^)} = . (52) 

The 5-function part has its origin from the virtual correction to the sphtting in the coUinear 
hmit. Pgg{z)\j^ is the same as our Pggi^z) defined in Eq. ()17|) except the 1/(1 — z) part has 
been regularized by the 'plus'-distribution prescription and so is replaced by 1/(1 — z)+. 
One can therefore write in a compact way the correction as 

Ct=o = -[ dzzPgg{z)\+ = [ dzPgg{z)\+ . (53) 



30 



2 _ 

This is a clever way to deduce the virtual correction from momentum conservation 
The last expression on the right is obtained by using the property Pgg{z) = Pgg{l — z) which 
follows also from the fact that in every branching two gluons are created for every one that 
gets destroyed. 

For gluon splitting in a medium because of the accompanying particle distributions, 
the virtual corrections cannot be deduced so simply especially when the system is out-of- 
equilibrium. However since the basic gluon splitting remains the same, one can reasonably 
expect that only the integrand of Eq. (j^Hj) will be modified by some factor involving the 
particle distributions of particle participating in the splitting. That is we expect 



1 



CT = --y^ dzPgg{z)WH{f,z), (54) 

where if is a function of the particle distributions and the momentum fraction z, which will 
be seen presently. 

Because of the medium as discussed in Sec. 1111 CI HTL corrections naturally establish 
limits on the z integration within which the coUinear approximation holds. These limits 
suffice as a regularization, we need not use the plus-distribution prescription. As a result 
the splitting function will always be written as Pgg{z) from now on. 



1. Type I: virtual correction with vacuum contribution 

We divide the virtual corrections into two types, one that has a contribution in the vacuum 
and the other has not. Bearing in mind that for every real process beyond the leading order, 
there is a virtual contribution that can be considered as its partner. They are related simply 
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FIG. 4: The self-energy insertion that gives the relevant medium virtual contributions. 

4 





FIG. 5: This figure shows the type I virtual correction on gluon 4 given by the expression in 
Eq. 



by originating from two different cuts of the same diagram |22|. For this reason the pair of 
real and virtual contribution are closely related especially in the coUinear limit. This facts 
allow them to be derived with some ease. The complete correction is a medium self-energy 
insertion shown in Fig. EJbut only the type I part of this will be discussed in this section. 

We start with the contribution that has a correction on gluon 4 of the 2-to-2 scattering. 
Very loosely this can be represented by 

9 'J 

X /l/2(l + /3)(1 + /4) - (1 + /l)(l + /2)/3/4]Fi(/4, • • •) • 

(55) 

Here -Fi(/4, . . .) to be given below carries additional particle distributions that arise because 
of the virtual correction and is a function of h and other internal variables of the second 
order a,np>.„de A.- Q. The a.pHtude car*, .he .ype I co„.ect,o. aUached 

to gluon 4. From the previous subsection we already know the form of the correction in 
the vacuum. Careful examination of the collinear limit of the virtual contribution shows 
that the structure surrounding the pole that gives rise to the collinear divergence in the 
vacuum is the same as the type I part of the medium contribution that gives rise to the 
same divergence. The only difference is whereas the vacuum part has a distribution of 1, 
the medium part has /. The only change is then in the appearance of a factor -F(/4, . . .) in 
the integrand of the integration over the phase space of the collinear gluons. The complete 
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contribution is 



C 



12^34 \Pl 



- j [C?P2]C?$34£'12,34^ ^1 '"^ 

+ /3)(1 + h) - (1 + + /2)/3/4] (l + f,{z) + U{1 



, (is asis) ^ , , 

12^34rX-_y dz-^P,,{z) 



(56) 

One can see in Eq. that there are two different factors that are made up of particle 
distributions. The first one in square brackets is the famihar difference of product of distri- 
butions that ensures this contribution to the collision terms will vanish in equilibrium. The 
second one in parenthesis [25 1 



Fi{U,z) = l + h{z) + h{l-z) 



(57) 



is a factor that comes from the virtual correction that is attached to gluon 4. They arise 
when an internal momentum loop has been converted into emission and absorption of gluon 
of the same momentum. For this type of correction, it is made up of the emission from and 
then followed by the absorption of gluons of the same momentum by gluon 4 [2^. This 
is depicted in Fig. El where the unlabelled lines represent the emitted and then absorbed 
gluons. The form of -F(/4, z) is deduced from detailed examination of this type I virtual 
correction in a medium which is done elsewhere j^]. Although no detail is given here, the 
correctness of this expression can be ascertained in the next section. 

In case that the virtual correction is on gluon 2, it follows similarly that 

,0 1 



^1(2)2 



12^34 



12^341 



X - / rfz^^P,,(.) 
2 / 3 2^^^^' 



If 1 X 

- / [rfp2]c?*34^12,34- \M 

/l/2(l + /3)(1 + /4) - (1 + /l)(l + /2)/3/4] (1 + f2{z) + ^(l - z) 



(58) 



2. Type II: virtual correction without vacuum contribution 



For the type II virtual corrections, these have no vacuum counterpart, when attached to 
jjluon 4, we have again loosely 



^ I [rfp2]t^$34/^12,34^ \ (-^12-34 ® + ^ 

.9 ■J 



12^34 



_^II(2)4 



^12->34 



X 



/l/2(l + /3)(1 + /4) - (1 + /l)(l + /2)/3/4 Fll(/4, • • •) • 



(59) 
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FIG. 6: This figure shows the type II virtual correction on gluon 2 and represents Eq. H62|) . On 
the left hand side, there is a spectator particle in the 12 — > 34 scattering. 

Here similar to the same in the previous subsection Fu is a function of and internal 
variables of the second order amplitude ^^12^34- Again the gluon splitting structure is 
the same independent of whether there is a medium or not, the only change when there 
is a medium is the factor of distribution Fji. After examining the structure of the virtual 
contribution, one can deduce 3] 

X /l/2(l + /3)(1 + h) - (1 + /l)(l + /2)/3/4] [hC-i^) - hC-)) . (60) 

The second factor of distribution in this case is 

Fu{Ux) = UC-^)-UC-) (61) 

which is a difference of two distributions because there is an absorption before an emission of 
a gluon of the same momentum. This is why this process can only occur in the medium. (See 
Fig.inifor a graphic representation when the correction is attached to gluon 2. The unlabelled 
lines on the right are the absorbed and then emitted gluons. The freely floating line on the 
left is a spectator particle for the interference to make sense.) The factor Fn(/4,x) is again 
somewhat mysterious but similar to the previous subsection, it comes from integrating out 
the internal loop from thermal field theory 3J]. This and the other results are deduced 
elsewhere [25 1. Even though we skipped a lot of steps to arrive at Eq. (jUUj) . one will see in 
the next section that the expressions in this subsection are correct. 

In case that the correction is attached to gluon 2, a similar expression can be derived in 
the same fashion 

-,11(2)2/ N / IJ-T. n I A ^ 1 2 



Cl2^34l/^l; 



12^341 



, dt asi\t\) „ , , 
t Zn 



— / [C?P2]C?$34^12,34:^ '"^ 
l^g J Z. 

fMl + h){l + h) - (1 + /l)(l + /2)/3/4] (/2(^) - /2(^ 



(62) 
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The graphs associated with this expression is shown in Fig. IHl 



IV. CANCELLATION OF THE COLLINEAR SINGULARITY BETWEEN THE 
REAL AND VIRTUAL CONTRIBUTIONS? 

We will now demonstrate that screening in an out-of-equilibrium medium is very impor- 
tant and more so than a plasma in equilibrium. This is because coUinear singularities cancel 
in the latter but not the former . Using the already derived results at the next-to- leading 
order, we split up the collision terms into the forward and backward process, for example 

Cl2^34: = C*12^34 — ^34^12 ■ (63) 

A. Both collinear gluons are in the final or initial state 

From Sec. IIII B II for real emission from final state particles, one gathers all the possi- 
bilities which are the three contributions: 3 || 4, 3 || 5, and 4 || 5. After relabelling they 
are 

C1II345 (Pi ) + C1II345 (Pi ) + ^£345 (Pi ) = 3 ^£345 (pi ) . (64) 

Explicitly this is 

X 1 f dz^^^P,,iz) 



^CStsM ^ Wrf$34/^12,34$^|Atl2^34p X^J 



S45 27r 



x/i/2(l + /3)(1 + Mz)){l + Ml - z)) . (65) 

In the vacuum (let's assume confinement has been turned off) one duplicates this same 
process by having gluon 1 colliding with gluon 2. The collinear divergence here must be 
cancelled by a virtual correction that has a vacuum contribution. Clearly the results from 
Sec. IIII D H are required. Summing the contributions with virtual correction on gluon 3 and 
gluon 4, after relabelling we have 

^1(2)3 ^1(2)4 _ ^ ^1(2)4 (f.ps 
'-'12^34 "T '-^12^34 — ^ '-"12^34 i \P^) 



where 



1 f ,_ d_sas{sl 

s 2n ' 



x/i/2(l + /3)(1 + U) (1 + U{z) + Uil - z)) . (67) 
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Comparing Eq. ()65|) with Eq. ()67|) . the only difference is in the factor of the particle distri- 
butions. That is 

Afip,) = (1 + U{z)){l + U{1 - z)) - (1 + + Uz) + Uil - z)} (68) 

(i) The plasma is in equilibrium: / takes the form of the Bose-Einstein distribution /be 
and A/^(p4) = is an identity, a consequence of the form of /be- This is the cancella- 
tion of collinear singularity that we mentioned in the introduction. This cancellation 
occurs even though these singularities are screened. A equilibrium system is doubly 
guarded against this type of singularity. 

(ii) The plasma is out of equilibrium: / can take on quite general forms and is only re- 
stricted by requirements such as the energy of the system must be finite etc. Therefore 
in general the above identity does not hold and A/^(p4) ^ 0. It follows that there 
is no longer any cancellation of the collinear singularity but since there is the failsafe 
mechanism of screening, there is no divergence even in the absence of any cancellation. 

In the case that both collinear gluons are in the initial state, the collinear absorption 
3 ^'315^12(^1) ^-^d the associated virtual correction 2 C^^l^^i2{Pi) to the 2-to-2 scattering again 
differ only by 

A/(P4) = f,{z)Ul -z)- /4{1 + /4(Z) + /4(1 - z)} . (69) 

After cancelling some of the terms in Eq. ()68p. this is exactly the same Ap{p4) factor above. 
The same discussion above clearly applies to collinear absorption too. 

B. One collinear gluon in each of the final and initial state 

In Sec. nil B 21 we saw that in the 12 345 5-gluon process there are the possibilities of 
collinear gluon pairs 2 || 3, 2 || 4 and 2 || 5. Collecting them together and relabelling, we 
have 

Cl Ji345 (Pl ) + '^lll345 (Pl ) + '^l1l345 (Pl ) = 3 C^l^^^^ {pi ) . (70) 

Explicitly this is 

X/l/2(i)(l + /3)(1 + /4)(1 + /2(i^)) . (71) 
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In the vacuum the colhnear divergence is cancelled after degenerate initial states are summed 
over. It does require other incoming gluons beside gluon 1 and 2 for that to happen. From 
the second type of virtual correction on gluon 2 to the 2-to-2 scattering in Sec. 1111 D 21 

x/l/2(l + /3)(l + /4)(/2(if^)-/2(i)) . (72) 

Comparing Eq. and Eq. (ff^ . the difference is entirely in the factor 

N'\P2) = /2(i)(l + f2{^)) - f2[h{^) - /2(i)] . (73) 

Once again we consider the effect of the state of the environment has on the interactions. 

(i) The plasma is in equilibrium: A/^^(p2) = is another identity when / is the Bose- 
Einstein distribution /be- This identity ensures the cancellation of any collinear loga- 
rithms between the real and virtual contributions. 

(ii) The plasma is out of equilibrium: the fact that the distribution can take many forms 
means that A/^^(p2) 7^ in general. However any hint of a divergence is rendered safe 
by having in the medium. 

Just as a gluon can be emitted from an initial gluon, so a final state gluon can also absorb 
an incoming one. Proceeding in a similar fashion as before, we find the difference between 
C!'^^^i2{Pi) C*34^'*^2(pi) is encapsulated in the factor 

A/"'(P2) = (1 + /2(i))/2(i^) - (1 + /2)[/2(i^) - /2(i)] = A/"(P2) • (74) 

This is the same expression as in Eq. (|7H|l and everything discussed above evidently also 
applies. 

V. IMPLICATIONS 

In the vacuum theory it is well-known that total cross-section is free of any collinear 
logarithms. This is true in an equilibrium system as well when / = /be as shown in the 
previous sections. The special form of /be ensures that the collinear logarithms are all 
cancelled. In neither of these cases have collinear processes any influence on the cross- 
section and the interaction rate is unaffected. However this paper is mainly concerned 
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with a system that is not yet thermahzed, as seen earher the colhnear logarithms remain 
even when the real and virtual processes are combined. To better appreciate what this 
means, let us consider the leading order 2-to-2 process together with only higher order 
collinear processes. In the vacuum and in an equilibrium system, the cross-section will only 
receive contribution from the 2-to-2 process in this special case. As soon as the system 
goes out of equilibrium, the cross-section will start receiving contributions at any order from 
the remaining collinear processes. There is a huge potential of significant changes in the 
cross-section or the gluon interaction rate in an out-of-equilibrium plasma especially when 
there is a large momentum transfer and that the transverse mass rrit is small enough so 
that approximately Q;sln(Q^/m^) ~ 1. Whether these collinear processes can help with 
thermalization is yet to be seen. However from the previous sections, it is clear that this 
depends largely on whether each pair of real and virtual collinear process summed up to 
give a positive or negative contribution to the collision terms. It was shown that the sum of 
each pair was controlled largely by the factors Af^{p) where a = I or II. The signs of these 
are in turned determined by the particle distributions. It is the state of the system that has 
the ultimate control: 

(i) it enables the collinear processes to contribute to the collision rate, 

(ii) it decides whether these collinear processes will help or hinder the thermalization of 
the system. 

The main question to be answered is if these collinear processes, which under normal 
circumstances would have no effect on the scattering rate, can have a helping hand in the 
equilibration of the system. From the previous section, the collinear contributions remain 
even after the sum of each real and virtual pair because of Af^i^p) ^ 0. Whether A/'^(p) > 
or Af"-{p) < will help determine largely how these non- vanishing contributions would have 
a positive or negative effect on the collisions (see below). If the equilibration of different 
momentum regions is completely random, then A/"(p) > in one region and Af°'{p) < 
in another there would be no significant net effect. Is there any reason for us to expect 
something different? 

In .f. .he e,„mb.a«on of a „.on pla..a . hea^ io„ wa. desCbed 

as a perturbative process powered by hard gluons. They function as an energy reservoir 
which pours energy into the softer momentum regions via gluon radiations so that the lower 
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energy gluons reach equilibrium first. This process gradually spreads upward from lower to 
higher momenta leading to the bottom-up thermalization scenario. This is reasonable for 
two reasons: (i) radiating multiple softer gluons is less restricted by energy conservation, 
(ii) the probability of emitting a lower energy gluon is bigger than emitting a higher energy 
one. As a result lower momentum regions that have an occupany level below the equilibrium 
level will be filled first ahead of any higher momentum regions. In this scenario one can 
expect the particle distributions at momenta below the range of the hard gluon reservoir to 
be organized in a way such that the particle occupancy at lower momenta is closer to the 
equilibrium value than at higher momenta up to just below the range of the reservoir of 
hard gluons. That is if we write 

Sfip) = f{p) - /be(p) , (75) 

then < —6f{p) < —6f{p') for any p < p', where p and p' are both below the range of the 
gluons that are powering the equilibration process. 

What impact has the bottom-up scenario on the results that we have shown so far? Be- 
fore we answer this question, let us remind ourselves what equilibration is. Obviously there 
are a number of ways that one can describe it, for example entropy production, momentum 
transfer between different momentum regions etc. Here we use the picture of particle occu- 
pancy in momentum space. Equilibration in this picture involves reducing the overoccupied 
momentum regions with respect to the equilibrium configuration by converting those gluons 
into gluons with momenta in the underpopulated areas. In heavy ion collisions one can 
think of the initial conditions as clusters of hard gluons that overpopulate hard momentum 
regions mostly aligned in the beam directions. Any other directions in momentum space 
would be underoccupied in this case. Alternately if the picture of the color glass condensates 
is used instead, the gluons would overpopulate the transverse plane in the central region and 
underpopulate anywhere else [3| • In both cases equilibration can be described as a process 
of transferring hard gluons via scattering and radiation from one direction or directions on 
a particular plane, which are overpopulated, to vastly different directions in other parts of 
momentum space that are underpopulated. 

Let us gather some important previously mentioned facts and prepare them for the central 
conclusion of this paper. These are the following: 
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(i) The possibility of equilibration via perturbative QCD is our starting premise for the 
simple reason that this has been well tested, better understood and also advocated in 



(ii) The most important contributions for equilibration are considered to be large an- 
gle collisions because they are efficient for transferring momenta between regions in 
momentum space. Yet by themselves they are not very probable so every hard scatter- 
ing should be augumented by small angle gluon radiations and/or absorptions. This 
augumentation is via the compensation of the small coupling by the large coUinear 
logarithm when ln((5^/'n^j) is of the order of unity. The large angle scattering thus 
acts to direct where the gluon radiations should go but it is the latter that actually 
fill out the underoccupied momentum regions. 

(iii) Because each of the collision terms is made up of the difference between the forward and 
the backward reactions, the most helpful processes in the collision terms for achieving 
equilibrium are those that can maximize this difference. There are two ways for this 
to happen. One is when all incoming gluons are from the overpopulated momentum 
regions and all outgoing ones are in the underoccupied ones. The other is having one 
incoming flux of gluons from the overpopulated region, the other incoming from an 
underoccupied region and all outgoing gluons are in the underoccupied ones. The 
underoccupied region that provides one of the incoming flux must be from a lower 
momentum region than the same from the overpopulated area. This is so that the 
higher occupancy of this lower momentum region can compensate partially for the 
larger backward reaction in this case. We will concentrate on scatterings that fulflU 
either condition. 

(iv) In the bottom-up equihbration scenario, the particle distributions in the underpopu- 
lated regions satisfy 




(76) 



for any p < p'. This implies that 



^f\p) = {i + f{m^ + fip-k))-{i + f{p)){i+f{k) + f{p-k)') 

= f{k)f{p -k)- f{p) (l + f{k) + f{p - A;)) 
> 0, 



(77) 
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and 

Af\p-k) = {i + f{k))f{p)-f{p-k)l^f{k)-f{p)') 

= (1 + f{p))f{k) - (1 + f{p - k)) (/(A;) - /(p)) 

= -{/(P - k)f{k) - f{p) (l + f{k) + f{p - k)) } 

= -A/^(p) 

< . (78) 

The simplest way to see that these relations are true is to consider a special case. We 
let f{k) = /be(A:), f{p -k) = /be(p - k) and f{p) = fMp) + ^fip) where 6f{p) < 0, 
so 

A/ \p) = -6f{p){l + /BE(fc) + /be(p - k)) = -Af\p -k)>0. (79) 

For more general distributions that satisfy Eq. (fTBj) . the relations Eqs. (f77j) and (|7Hj) 
can be more readily verified numerically. 

(v) In the overpopulated region where the hard gluons provide the energy for the equili- 
bration, the lower energy gluons here would be closer to equilibrium than the higher 
energy ones so we expect for most of these hard gluons 

Sfip) > Sfip) > (80) 

for p < p' . At the tail of the distribution at very high energies, this is no longer 
true but the occupancy here becomes too low to have too much significance. The 
consequence is that the relations of the A/'s in (iv) are reversed in the hard gluon 
reservoir 

M\v-k) = -M\v)>^- (81) 



With these points in mind, we now apply them to the gluon interactions. Without loss 
of generality, we let p\ to be in the hard gluon reservoir which supplies the energy for the 
equilibration. Up to the next-to-leading order, the collision terms with the inclusion of only 
collinear processes at and up to order are 

C(pi) = Ci2^34(pi) + {iCt^^M + 2CSll54(pi)) + (4Cf|U5(pi) + 2CS4(Pi)) 

+2C^ilL45(Pi) + 3C||l345(pi) + C'lS^M + C^£34(Pi) + 0{al) . (82) 
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Those terms in parenthesis are the pairs of real and virtual contributions. Explicitly they 
are 



2^1(2)4 



12^34 (Pl 



j Mrf$34^12,34 5^|A^12^34|' X dz^^^Pgg{z) 



X 



/l/2(l + /3)-(l + /l)(l + /2)/3 Afip,) 



(83) 



-— / [rfp2]c?*34^12,34 J^l^ 



12^341 



, dt as(\t\) „ , , 



t 27r 



X 



/l/2(l + /3) - (1 + /l)(l + /2)/3 A/"(p4 



(84) 



and 



3C 



12^345 (Pl) +^{£34(^1) 
— / [C/P2]C^$34^12,34:T ^ \M 

+ /3)(1 + U) - (1 + /l)/3/4l A/"(P2 



12^34 



9 ; , dt as{\t\) ^ , , 



C 



2||3 

123^45 



- / [rfp2]rf$34^12,34- 1-^ 
/l(l + /3)(1 + /4) - (1 + h)hf^N\V2) . 



12^34 X 



s 27r 



(85) 



(86) 



Note that each of these expressions carries either A/^ or A/^^ as the last factor. These 
expressions are possible because A/° have at least two equivalent forms. One is associated 
with the forward reactions and the other with the backward ones but they are in fact 
equivalent after some algebra as seen in Eqs. (|77jl and (fTHjl . As they stand these expressions 
are hard to analyze, we elect to consider the scattering cross-section instead. As a result 
only the forward part of these equations will be used below when we try to construct the 
scattering cross-section. 



A. The in-medium cross-section 



We consider the interactions in the medium initiated by two flux of incoming particles. 
One flux consists of coUinear particles centered on the momentum pi. By this we mean gluons 

29 



of momenta pi, or collinear gluons with momenta that can sum up to pi such as pi — p and 
p where pi \\ p, |pi| > |p| etc are included in the flux. The second flux is centered on p2. 
Through the study of the collisions produced by these two flux, we can deduce the effect of 
the collinear processes on the collision rate or equivalently the interaction cross-section 

.>„..W..(P,.P.) = -M!(2p»%M + 2rf%M)^ (87) 

Ug V d-^p2 d-^pi } T\T2 Vr 

where J^j, i = 1,2 are the particle flux and Vr is the relative velocity between the two 
flux. The expression of C^^jr^ denotes those forward reactions in the collision terms that 
have contributions from interactions between the pi-centered with the p2-centered flux and 
the associated virtual corrections. This definition of the cross-section reduces to the usual 
expression when in the vacuum. The explicit expression C^^j:^ up to o?^ can be extracted 
from Eq. ©, Eqs. (|HS1), dHHl), (Ell) and ^ 

C^.^M) = ICu^uiPi) + \ (3C£345(Pl) + 2CS4(Pi)) 

Z \ / \ / modified 

+ (cS_45(Pi) + C^SlPi)) • (88) 

The factor of one-half in the first three terms is to avoid including those interactions that are 
basically symmetric in gluon 1 and 2 twice in a"in-medium- C^_^jr^{p2) has basically the same 
expression except gluon 1 is interchanged with gluon 2. The two expressions C^_^jr^{pi) and 
C^^jr^{p2) together cover all the interactions in C{pi) except C^-^-*45ip^) ■ This particular 
interaction is initiated by gluon 2, 3 and not 1, thus it does not contribute to crjn- medium (Pi,P2)- 
Note that the expression Eq. (j85|) have to be modified first because the incoming momentum 
in this expression is roughly P2/X instead of p2- After some adjustments including relabelling 
V2 v'2 — ^^2, gluon 3 and gluon 4 to distinguish them from that in Eqs. (|H^. (jHU)) and 
and freezing the value of the incoming p2 instead of if becomes 



I modified 



3^^lE345(Pl)+C^S34(Pl) 
^ y"Mt^$3'4'/^12',3'4'^$]|-Ml2'^3'4f X ^ t/xy ^^P,, (x) 

+ /30(1 + /4') - (1 + /l)/3'/4'l N\W2) . (89) 



The d/d^Pi in Eq. (jH7|) is there to indicate that the integration over d^pi in the different C's 
should be removed and the 4-momentum at Pi is frozen at that value algrebraically. 
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Let us now examine each term beyond the leading order term in Eq. (j88j). Here a term 
at al refers to each pair of real and virtual contribution enclosed in parenthesis. From the 
above discussions, point (iv) guarantees that the second and third term give a positive and 
negative contribution respectively to o"in-medium- These two terms differ only by a collinear 
emission 4 — ij from and a collinear absorption 4j i hj one of the final particles. Apart 
from the ds and dt integrations that yield the slowly varying logarithms Eqs. ()38p and ()45|) 
which are both essentially controlled by the Q'^/m'f ratio, one has the dz integration over 
|A/^(p4) and the other has dx over A/^^(p4). While A/^(p4) > — A/^^(p4) > for all in 
the underpopulated momentum region, the same cannot be said of |A/^(p4) > — A/"(p4) 
which is not always true. However 

lim A/ \p^) = hm A/ ^(^4) — > - lim A/"(p4) • (90) 

While Ap^lpi) decreases monotonically towards zero as x decreases from unity, A/^(p4) has 
a minimum at 2; = 1/2 and two maxima at z = and 2; = 1 as can be seen in Fig. [3 These 
maxima can be quite large. This leads to the result that after integrating over the light- 
cone fractions x and z in the second and third terms of Eq. ()88|). the second and positive 
contributing term will be larger than the third negative term (see Fig. [Tj). The sum of the 
second and third terms then contribute positively to CTin-medium- Of course the limits that 
should be used are z^in < z < z^g,^ and Xmin < x < Xmax shown in Sec. IIII (^I but these are 
not too far from and 1 when the ratio of Q'^/mf is large. 

There are the remaining fourth and fifth terms. For these we have to distinguish the two 
different incoming gluons raised in point (iii). 

1) Both flux from overpopulated momentum regions: 

The fourth term of Eq. §^ has Af\xp2) = -Af\p2) > from Eq. ^ (point (v)), 
as part of the integrand and the fifth has Af^{p2) from Eq. ()86|) . Thus the fourth has a 
positive contribution to the cross-section and the fifth contributes negatively. Although 
they both carry the same factor A/^(p2), the fifth term has an extra symmetry factor 
of one-half and the occupancies of the outgoing particles are smaller than those from 
the fourth term because the outgoing gluon momenta are in general larger. To see 
this, the outgoing momenta of the fourth term are p^i + ^4/ = Pi + P2 — Pb whereas 
those of the fifth term are Pz + Pi = Pi+P2- f{Pw) and f{pii) will in general be larger 
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FIG. 7: A comparison of the modulus of the typical Af (p) and A/ (p) for some fixed p in an 
underpopulated momentum region for various values of the lightcone fraction z or x. The area 
covered by the curve of A/^(p) is clearly more than twice that of A/^^(p). 

than /(ps) and /(P4) as a consequence. The sum of the fourth and fifth term therefore 
gives a net positive contribution to cTiji.jncdium- Grouping now all the contributions of 
coUinear processes from the final states as well as those in the initial states discussed 
here, the sum of the four next-to-leading contributions in Eq. (jHH|) are negative so their 
contribution to the in-medium cross-section is positive. 

2) One flux from an overpopulated and the other from an underpopulated region: 

In this case the p2-centered flux is from an underpopulated region since pi has already 
been chosen to be from the hard overpopulated region. This is the reverse of situation 
of 1). From point (iv) the fourth term of Eq. (jHHj) is now negative with — A/^(p2) < 
and the fifth becomes positive with /\f^{p2) > 0. The same reasons given in 1) that 
ensured the sum of the two was positive there now ensure that it is negative. To sum 
it up, the second and third term of Eq. (jHH|l give a net positive contribution but the 
fourth and fifth term give a net negative contribution. 

Of course these are not the only contributions because both C^_^jr^{pi) and C^-^jr^iP2) 
have an entry in (Tin-medium in Eq. (IHTjl . In the Cp^jr^{p2) entry, the second and third 
term and also the fourth and fifth term as well would still give net positive contributions 
because there is not much difference here from case 1). It follows that the coUinear 
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processes in the initial states of the two incoming flux in this case tend to have opposing 
and compensating effects on each other. In general it is hard to tell which contributions 
will have the upper-hand since it depends to a certain extent both on the value of pi 
and p2, and on how close each of these is to the boundary separating the overpopulated 
to the under-occupied regions. Thus the coUinear processes in the initial states do not 
have much net effect on the average, however those in the flnal states are still very 
much present and positive. One can conclude that there is still enhancing effect on 
c"in-medium but to a Icsscr degree than that in 1). 

VI. IN CONCLUSION 

Because of the overall positive contributions from the coUinear processes together with 
the fact that these contributions vanish when the interactions occur in the vacuum or when 
the system is in equilibrium, it means that the in-medium non-equilibrium cross-section is 
larger or equivalently the out-of-equilibrium interaction rate is higher than what is expected. 
This is the novel perturbative mechanism for achieving very fast thermalization within 0.6 
fm/c at RHIC that we are hoping for. Underlying this mechanism is the non- vanishing of 
coUinear logarithmic terms in a system that is out of equilibrium. We are not aware that 
this has been pointed out in any existing literature. This is considered in more detail in 
[2^. In closing we have yet to consider here the size of the enhancement which is clearly 
dependent on how far away from equilibrium the system is, and also we have not considered 
how the contributions at higher orders will affect our main conclusion, nevertheless a lot of 
the structures in the various C(pj)'s, the factors of distributions and the inequalities that 
they satisfy are shared with those at higher orders. These are better studied and verifled 
numerically than by using the arguments that we have attempted here. All these will be 
done in the future. 
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